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The anomalous X-ray pulsar 4U 0142+61 was observed with Suzaku on 2007 
August 15 for a net exposure of ~100 ks, and was detected in a 0.4 to ~ 70 keV energy 
band. The intrinsic pulse period was determined as 8.68878 ±0.00005 s, in agreement 
with an extrapolation from previous measurements. The broadband Suzaku spectra 
enabled a first simultaneous and accurate measurement of the soft and hard com- 
ponents of this object by a single satellite. The former can be reproduced by two 
blackbodies, or slightly better by a resonant cyclotron scattering model. The hard 
component can be approximated by a power-law of photon index Fh ~ 0.9 when the 
soft component is represented by the resonant cyclotron scattering model, and its 
high-energy cutoff is constrained as > 180 keV. Assuming an isotropic emission at 
a distance of 3.6 kpc, the unabsorbed 1-10 keV and 10-70 keV luminosities of the 
soft and hard components are calculated as 2.8 x 10^^ erg s~^ and 6.8 x 10^^ erg s~^, 
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respectively. Their sum becomes ~ 10'^ times as large as the estimated spin-down 
luminosity. On a time scale of 30 ks, the hard component exhibited evidence of 
variations either in its normalization or pulse shape. 

Key words: magnetic fields — neutron stars: individual (4U 0142+61) — X-rays: 
general, individual (4U 0142+61) 

1. Introduction 

Anomalous X-ray pulsars (AXPs), comprising at present ~9 objects discovered in the 
local universe, are characterized by rotational periods in the range P = 5 — 12 s, and spin down 
rates as P ~ 10"^^ s s~^. Together with soft gamma repeaters (SGRs), AXPs are thought 
to form a subgroup of neutron stars called "magnetars", of which the surface magnetic field 
strengths are believed to reach ~ 10^^~^^ G (Thompson & Duncan 1995; Thompson & Duncan 
1996; Woods & Thompson 2006). In energies below ~ 10 keV, AXPs exhibit spectra which are 
much softer than those of ordinary binary X-ray pulsars, while harder than those of isolated 
neutron stars. Their X-ray luminosities (10^^~^^ erg s~^), which often exceed up to two order 
of magnitude those available from their spin-down power (~10^^ erg s~^), are usually thought 
to be sustained by a release of energies stored in their ultra-strong magnetic fields (Duncan & 
Thompson 1992; Thompson & Duncan 1995). 

A novel observational window onto magnetars has been opened by the INTEGRAL 
discoveries of a distinct pulsed hard X-ray component from at least 3 AXPs and 2 SGRs 
(Kuiper et al. 2004; Mereghetti et al. 2005; den Hartog et al. 2006; Kuiper et al. 2006; Gotz 
et al. 2007). This component extends to ~ 100 keV or more with a very flat photon index of 
r ~ 1, and exhibits a luminosity comparable to or higher than that in the softer X-ray band. 
Although theoretical accounts for this enigmatic component are far from settled (e.g., Heyl & 
Hernquist 2005; Beloborodov & Thompson 2007; Baring & Harding 2007), its near-absence in 
other types of X-ray sources suggests its close relation to the proposed strong magnetic fields 
of magnetars. 

The study of magnetars now requires a broad energy band, because their soft and 
hard components, which are generally variable, must be measured simultaneously and accu- 
rately. This makes Suzaku (Mitsuda et al. 2007) an ideal observatory. In fact, the X-ray 
Imaging Spectrometer (XIS), operating in the 0.2-12 keV range (Koyama et al. 2007), provides 
high-quality data of their soft components. Simultaneously, the Hard X-ray Detector (HXD; 
Takahashi et al. 2007; Kokubun et al. 2007), consisting of HXD-PIN (10-70 keV) and HXD- 
GSO (40-600 keV), can conduct detailed spectroscopic studies of their hard components in 
considerably shorter exposures than are needed by INTEGRAL. 

Such wide-band spectroscopic observations with Suzaku have allowed detections of the 
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hard-tail component from activated magnetars, including SGR 0501+4516 (Enoto et al. 2009; 
Rea et al. 2009; Enoto et al. 2010b) and IE 1547.0-5408 (Enoto et al. 2010a). Through these 
observations, the hard X-ray emission from SGR 0501+4516 was revealed to form a component 
distinct from the soft X-rays, while the power-law-shaped hard-tail component of IE 1547.0- 
5408 to extend to the soft X-ray band below 10 keV. Thus, a key point of the observation of 
magnetars is to detect and quantify their wide-band spectra simultaneously, without hampered 
by cross-calibration uncertainties or non-simultaneity involved in multi-satellite observations. 
This is a great advantage provided with Suzaku. 

In the present paper, we report on a high sensitivity broadband observation of a pro- 
totypical magnetar 4U 0142+61 made with Suzaku. This object is one of the most luminous 
AXPs with P~8.68 s (Israel et al. 1994), located at the Galactic anti-center direction, probably 
on the Perseus arm (Durant & van Kerkwijk 2006a). Its pulsed hard X-ray emission was 
detected up to ~ 200 keV by INTEGRAL and RXTE (Kuiper et al. 2006; den Hartog et al. 
2006; den Hartog et al. 2008). On the other hand, past COMPTEL observations gave a flux 
upper limit in energies above 750 keV (Kuiper et al. 2006; den Hartog et al. 2006). If the source 
exhibits little time variations between these observations, the hard X-ray spectrum must turn 
over at energies between ~200 and ~750 keV, though it has not yet been clearly detected. 

2. Observation and Data Reduction 

We observed 4U 0142+61 with Suzaku from 04:04 UT on 2007 August 13 to 12:30 UT 
on August 15, for a gross duration of 203 ks and a net exposure of ~ 100 ks. The XIS and the 
HXD were both operated in their normal modes, except that 1/4 window option was applied 
to the XIS to achieve a time resolution of 2 s. 

We placed the source at the HXD nominal position, and adjusted the roll angle to avoid 
contaminating bright hard X-ray sources. In particular, we tried to minimize the confusion 
with the Be X-ray binary RX J0146.9+6121. Located 24' off 4U 0142+61, this source exhibited 
a 20-50 keV flux comparable to that of 4U 0142+61 during INTEGRAL observations in 2003- 
2004 (den Hartog et al. 2006). As a result of the roll-angle adjustment, the HXD-PIN effective 
area onto RX J0146.9+6121 was <29 % of that onto 4U 0142+61 in the 10-70 keV energy band. 
Based on this, and the reported steep spectrum (F ~3; den Hartog et al. 2006), we regard the 
contamination of RX J0146. 9+6121 to the HXD negligible. There were no other catalogued 
bright X-ray sources inside the HXD field of view. 

We used the Suzaku data prepared via version 2.1 processing, and the HEASOFT version 
6.4 or later tools. The employed data screening criteria were; (a) the time after an exit from the 
South Atlantic Anomaly should be >500 s and >436 s for the HXD and the XIS, respectively, 
while the time to the next entry should be >180 s for the HXD. (b) the target should be above 
the Earth rim by >5° for the HXD and the XIS, and, in case of the sunlit Earth rim, that 
should be >20° for the XIS; (c) the instantaneous pointing direction should be within 1'.5 of the 
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mean; and (d) the geomagnetic cutoff rigidity should be >6 GV. As a result of these screenings, 
we archived a net exposure of 99.7 ks with the XIS, and 94.7 ks with the HXD. 

3. Data Analysis and Results 

3.1. Background subtraction 

We accumulated the screened data of the XIS from a region within 2'.0 (2.8 mm) radius of 
the source centroid, and derived a background spectrum from another region which is symmetric 
with respect to the XIS field-of-view center. Corresponding rmf and arf files were created using 
xisrmfgen and xissimarf gen tools (Ishisaki et al. 2007), respectively. 

The non-Xray background (NXB) of the HXD must be subtracted using appropriate 
models that emulate NXB events to be observed during the on-source exposure. Specifically, 
we utilized LCFITDT (bgd_d 2.0 ver0804) and LCFIT (bgd_d 2.0 ver0804) methods (Fukazawa 
et al. 2009) to create the fake NXB events of HXD-PIN and HXD-GSO, respectively. These 
fake photons were also screened by the same criteria as for the on-source event data. 

Using the data acquired during Earth-occultation periods (~20 ks in total) in the present 
observation, we confirmed that the HXD-PIN and HXD-GSO spectra obtained therein can be 
reproduced by the employed NXB models both within which agrees with the typical 

uncertainties of these models. This calibration suggests that we are slightly (~0.6%) under- 
subtracting the NXB of HXD-GSO. Thus, we finally subtracted the HXD-GSO background 
after increasing it by 0.6%, and include systematic uncertainties of 0.6% after Fukazawa et al. 
(2009). 

3.2. Source detection 

Figure 1 shows background-subtracted light curves of XIS3 and HXD-PIN. For the soft 
X-ray light curve, we utilized the XIS3 data because those from XISO or XISl are more subject 
to artificial intensity variations synchronized with the orbital period of Suzaku, caused by 
thermal distortions of the spacecraft structure (Uchiyama et al. 2008). We do not show a light 
curve of HXD-GSO, since the NXB model for HXD-GSO is valid only for integrations longer 
than ~100 ks for such low count rate objects. Thus, the source has been detected clearly 
not only by XIS3 but also by HXD-PIN, as the average HXD-PIN count rate, ~ 0.08 c s'^ 
significantly exceeds the expected CXB rate of ~ 0.03 c s~^ (Kokubun et al. 2007). 

At the binning as shown in Figure 1, the XIS3 count rate was constant during the 
observation within 15%. The possible variations seen in figure la can be mostly attributed to 
slight losses of signal photons, due to the thermal distortion of the spacecraft structure which 
caused the image outskirts to partially fall outside the 1/4 window. The HXD-PIN signals, 
with 550 s binning, are also roughly constant: the NXB-subtracted light curve can be fitted 
by a constant with xV^ = 154.9/126. However, the HXD-PIN light curve may be slightly 
decreasing with time, because it can be fitted better (j^ jv = 142.2/125) by a linear function 
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where the count rate at the end of the observation is 35% lower than that at the beginning. 
Both the XIS3 (2-s bin) and HXD-PIN (1-s bin) counts are consistent with obeying Poisson 
distributions, without any significant short-term increases of their count rates. 

Figure 2a shows background-subtracted spectra of 4U 0142+61 obtained with three XIS 
cameras, HXD-PIN, and HXD-GSO. The background levels are also shown. From the HXD- 
PIN spectrum, we subtracted not only the modeled NXB described in section 2, but also the 
CXB spectrum taken from Moretti et al. (2009). Although the CXB brightness within the 
HXD-PIN field of view is expected to vary by ~ ±12% from sky to sky (Fukazawa et al. 2009), 
the effect is of the order of 0.5% of the NXB, and hence is negligible compared to the NXB 
modeling uncertainty. The CXB contribution itself is negligible in the HXD-GSO data. 

As seen in figure 2a, the HXD-PIN signals (after removing the NXB and CXB) are 
statistically significant over the entire 10-70 keV band, and stay > 10% of the NXB level 
therein. Since the NXB can be reproduced typically down to ^2% (Fukazawa et al. 2009), 
the HXD-PIN detection is highly significant with ^3ct even considering the systematic errors. 
In fact, this NXB modeling uncertainty is considerably smaller than the statistical errors. In 
contrast, the NXB-subtracted HXD-GSO signals are found at a level of ~ 1% of the NXB 
which is a typical systematic uncertainty. Therefore, the source detection with HXD-GSO is 
considered insignificant only from the spectrum. 

In order to grasp basic properties of the spectra of 4U 0142+61 in a model independent 
way, we normalized them to those of the Crab Nebula (a power-law with a photon index P = 2.1). 
We utilize the HXD data of the Crab Nebula obtained on 2007 March 20, while employ the 
simulated Crab Nebula spectrum for the XIS, in order to avoid pile- up effects at the source 
centroid. The derived "Crab ratio" is presented in figure lb. Thus, in addition to the long- 
known soft component, the hard component (section 1) is very clearly visible, rising up from 
~10 keV and reaching 10 mCrab at 100 keV. 

3.3. Source pulsations 

After applying barycentric corrections to the arrival times of individual photons (Terada 
et al. 2008), we searched the background-inclusive XIS and HXD data for the 8.7-s source 
pulsation via standard periodogram analyses. Since the spacecraft wobbling is insignificant on 
the relevant time scales, the data from XISO and XISl ware also incorporated. 

Figure 3 show periodograms obtained with the XIS and HXD-PIN in the 0.8-10 keV 
and 12-40 keV energy range, respectively. Thus, the pulsation has been detected clearly by the 
XIS, at a barycentric period of 

P = 8.68878 + 0.00005 s. (1) 

This agrees, within 1 1 fis, with a value extrapolated from previously measured P and the period 
change of P = 1.960 x 10"^^ s s~^ (Gavriil & Kaspi 2002). The width of the XIS periodgram peak 
is ~ 3 X 10~^ s, or ~ 4 X lO^^P. This is consistent with the gross time span of our observation, 
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~ 2.0 X 10^ s (section 2), wherein ~ 3 x 10^ pulses are contained. 

In the 12-40 keV HXD-PIN periodgram, a peak with x^l^ = 35.9/15 = 2.4 is found at 
a consistent period. This gives a chance probabihty of 8.4 x 10~^ for the HXD-PIN profile 
to be caused by Poisson fluctuations. An independent and bin-free Z^-test (Buccheri et al. 
1993; Brazier 1994), together with a widely used choice of = 2, gives = 21.8 for the HXD 
data at the period of (1), with a chance probability of 2.1 x 10~^. These tests reveal that the 
HXD-PIN signals are pulsed, at the period of equation 1, with a high significance (>99.9%), 
even though this level would not be high enough to support the detection of a periodicity 
without the aid of the XIS data. 

Fgure 4 shows pulse profiles in the 0.8-4.0, 4.0-10, 10-70, and 80-150 keV energy ranges 
observed the XIS (XISO, XISl, and XIS3), HXD-PIN, and HXD-GSO data, all folded at the 
period of equation (1), employing an epoch of 54326.356 (MJD). Thus, the pulsation is less 
significant in the 80-150 keV HXD-GSO data, with xV^ = 12.2/6 = 2.0 (null hypothesis proba- 
bility 0.06). The pulse profiles exhibit some energy- dependent changes, in agreement with pre- 
vious reports (Kuiper et al. 2006; den Hartog et al. 2008). The peak-to-bottom pulsed fraction, 
(-^max — -Fmm)/(-^max + -^min), where Fmax and Fmin are the maximum and minimum background- 
subtracted count rates across the pulse phase, are evaluated to be 8.4±0.3%, 13±1%, and ~41% 
in the 0.8-4, 4-10, and 10-70 keV energy ranges, respectively. 

3.4- Short-Term Variations 

Further periodogram analyses of the HXD-PIN data gave an indication that the pulsa- 
tions in the 12-40 keV band become more prominent toward the latter part of the observation. 
We hence divided all the exposure into three epochs of ~30 ks each, and produced three sets of 
12-40 keV HXD-PIN light curves folded at the period of equation (1). As shown in figure 5, the 
pulse profiles from the 1st and 2nd epochs were found to be consistent with each other within 
statistical errors. In contrast, as shown in panels (a) and (c) of figure 5, those from the 1st and 
last epochs appear considerably different; the pulses became deeper and more double-peaked in 
the 3rd epoch, while the average intensity decreased in agreement with figure 1.1 These effects 
are unlikely to be caused artificially by cumulative phase errors, which is estimated to be <9% 
of a cycle when folded at equation (1). 

To quantitatively examine the suggested pulse-profile changes, we took phase-by-phase 
differences between the 1st and 3rd epoch pulse profiles, and obtained = 20.8 for u = 9. The 
calculation took into account (in quadrature) the typical systematic background uncertainty 
by 0.004 c s~^ (1 sigma). Since the implied null- hypothesis probability is 1.3%, the two folded 
profiles are inferred to be different at a 98% confidence level. This suggests that the 12-40 keV 
signals changed, either in the average intensity or in the pulse profiles (or both). However, these 
alternative possibilities cannot be distinguished, since the significance of the profile difference 
reduces to = 15.3 for i/ = 9 (or a chance probability of 8%) if the difference profile is fitted 



6 



by a constant. 

For reference, the XIS pulse profiles in 0.4-4 keV did not change by more than 7% among 
the 3 epochs. Any spectral change in the XIS or HXD-PIN band was not significant among the 
3 epochs, either. 

3.5. Spectral Analysis 

3.5.1. Wide-band spectra 

To quantify the spectral soft and hard components simultaneously, we jointly fitted the 
data from the XIS, HXD-PIN, and HXD-GSO, in the 0.8-10, 12-70 keV, and 50-200 keV ranges, 
respectively. Although the signal detection with HXD-GSO is not significant, we include the 
GSO data points with the systematic error properly added, because we can use them at least 
as upper-limit points. The on-source and corresponding background spectra were prepared as 
stated in section 3.1. We also incorporated the correction factor for HXD-GSO spectra, which 
has been introduced to reproduce the Crab spectra by a single power-law (Takahashi et al. 
2008). 

As figure 2b shows, the XIS and HXD data points happen to define, in a rather clear-cut 
way, the soft and hard components, respectively. As the simplest attempt, we separately fitted 
the XIS spectra with an absorbed blackbody model and the HXD one by a power-law model, 
to find that the latter is successful with Fh = 0.8lto;2 (the former is not; see below). As a next 
attempt, we tried to jointly reproduce the whole data simultaneously by an absorbed blackbody 
model plus a power-law model, representing the soft and hard components, respectively. The 
model normalization between the XIS and HXD-PIN was constrained to be 1.18 (Maeda et al. 
2008). When the hard power-law photon index is fixed, for simplicity, at the above value of 
Fh = 0.8, this model gave a blackbody temperature of kT ~ 0.46 keV. However, as shown in 
figure 6c, the fit was unsuccessful with = 7326.9/412 = 17.8, due to significant residuals 
from the blackbody model at ~1.2 keV and ~5 keV. The resultant parameters are shown in 
table 1 as Model A. When F is allowed to be vary, it became very steep (F ~ 3.7) so as to 
account for the model inaccuracy below ~10 keV rather than to reproduce the hard tail. 

Trying to improve the fit, we added a second steep power-law, so that the soft component 
is now expressed by a sum of a blackbody and a power- law (i.e., BB-I-PL). The obtained best-fit 
solution is described by a blackbody of kT = 0.43 keV, a steep power-law of F = 4.0, and a hard 
power-law of Fh = 0.1 (which is now left free), with a three times larger A'^h value than Model 
A. Systematic errors of F^ are estimated to be ~30% when considering the 2% uncertainty 
of the PIN-NXB. However, as shown in figure 6d (Model B in table 1), this model was still 
unsuccessful with x^/^ = 589.9/409 = 1.44, mainly due to deviations in the low energy range 
of the XIS. 

By replacing the steep power-law with a second blackbody, and thus representing the 
soft component by two blackbodies (or 2BB), the fit was significantly improved to x"^/^ = 
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459.8/409 = 1.12; these results are shown in figure 6e (Model C in table 1). The two blackbody 
temperatures were obtained as kTi = 0.34 keV and kT2 = 0.63 keV, with their emission radii 
(assuming circular regions) as 13.2 km and 2.3 km, respectively, when assuming a distance 
to the source as 3.6 kpc (Durant & van Kerkwijk 2006a). The absorbing column became 
Nn = 6.4 X 10^1 cm-2. 

Although the above Model C is approximately successful, the fit still leaves significant 
residuals. As one of alternative empirical models for the soft component, we employed, after 
Tiengo et al. (2005) and Enoto et al. (2010c), the Comptonized blackbody model described in 
Enoto et al. (2010a). As shown in figure 6e and summarized as Model D in table 1, this model 
gave a better fit; x'^/'^ = 441.5/410 = 1.08, with kT = 0.29 keV, a soft power-law photon index 
Tg = 3.8, and rh=0.4. The results using this particular model was already reported previously 
(Enoto et al. 2010c). 

In the strong magnetic field as ~10^^ G, resonant cyclotron scattering of soft X-ray 
photons becomes more efficient than Thomson scattering. This effect can be expressed by 
"Resonant Cyclotron Scattering (RCS)" model (Thompson et al. 2002; Lyutikov & Gavriil 
2006; Rea et al. 2008), which calculates a blackbody spectrum modified by resonant cyclotron 
up-scatterings off hot electrons in the neutron star magnetosphere. This model, combined with 
the hard power-law, has given a still better account of the data over the 0.8-200 keV energy 
band, with x'^/'^ = 362.0/409 = 0.88. This fit is shown in figure 6a and figure 6f, and in table 1 
as Model E. The derived blackbody temperature (~0.3 keV) is similar to that of Model E, while 
the Compton optical depth, which was left implicit in Model D, became r ~ 2.0. The photon 
index of the hard power-law became Fh = 0.89lQjQ(stat.)^o;}o(sys.), which is close to the value 
of Eh ~ 0.80 obtained using the HXD-PIN data only. 

Figure 7 represents the uFi, spectra of 4U 0142+61, fitted with the most successful RCS 
plus power-law model (Model E) described above. The X-ray fiux in the 10-70 keV energy band 
is 4.4 X 10~^^ erg cm~^ s~^, while the unabsorbed 1-10 keV fiux is 1.8 x lO"^'' erg cm~^ s~^. If 
the hard component extends to 200 keV, its fiux in the 10-200 keV band becomes 1.5 x 10^^" 
erg cm~^ s~^, which is comparable to that of the soft component. 

3.5.2. Possible local spectral features 

At a closer inspection, we still find negative residuals at ~ 20 keV and ~ 40 keV in 
the HXD-PIN spectrum (figure 2). Since the NXB spectrum of HXD-PIN shows no spectral 
features at ~ 20 keV, and only a very weak Gd-K„ emission line at ~ 42 keV, the suggested 
features cannot be instrumental. However, we refrain from further analysis of these features, 
because their statistical significance is not high enough. 

Lastly, we examined the hard power-law for a high-energy cutoff, because the hard-tail 
component must steepen at some energies, in order for its luminosity not to diverge. For this 
purpose, we employed Model D, and multiplied its hard power-law with an exponential cutoff 
factor of the form exp(— E'/E'cut), where E is the photon energy and E^ut is a cutoff parameter. 
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Then, the data gave a lower hmit as Ecut > 180 keV at the 90% confidence hmit. The choice of 
Model D is because it gives (among the successful 3 models) the fiattest value of Fh, and hence 
the derived lower limit on Ecnt is considered to be most conservative. 

4. Discussion 

4.1. Wide-band spectra of magnetars 

Recent X-ray observations of magnetars show that their persistent broad-band spectra 
are commonly represented by an optically thick soft component of kT ~ 0.3 keV, and a power- 
law-like hard-tail component which dominates in energies above ~10 keV. Besides, it was 
revealed that the two-component spectrum depends significantly on their characteristic age, 
in such a way that the hard-tail component becomes weaker (relative to the soft component), 
yet harder, towards sources with older characteristic ages (Kaspi & Boydstun 2010; Enoto et 
al. 2010c). When systematically studying such spectral properties, including in particular the 
evolutionary behavior revelead by Enoto et al. (2010c), 4U 0142+61 is considered to play a 
particular important role, for the following three reasons. 

First, among the magnetars in the evolutionary study by Enoto et al. (2010c), 
4U 0142+61 has the oldest characteristic age, 70 kyr, except for IE 2259+685 of which the 
hard-tail component is not well constrained. As a result, 4U 0142+61 is currently the old- 
est magnetar with the well determined hard-tail component, and hence it significantly rein- 
forces the age-correlated spectral evolution of magnetars (Enoto et al. 2010c). Secondly, as a 
common characteristic of aged magnetars (Enoto et al. 2010c), the two spectral components 
of 4U 0142+61 are well separated from each other (see figure 2). This, combined with the 
relatively high X-ray intensity of 4U 0142+61, allows us to better quantify its two spectral 
components than in other objects of this class. Finally, the extreme fiatness (Fh ~ 0.8) of the 
hard tail of 4U 0142+61, which is again typical of aged magnetars, is expected to provide very 
strong constraints on the emission mechanism of the magnetar hard X-rays. 

4-2. Soft X-ray emission of 4U 0142-1-61 

When comparing among the soft X-ray spectral models in § 3.5.1, the blackbody plus soft 
power-law model (Model B) requires a higher column density of photo-absorption than the other 
models. This is likely to be an artifact, necessitated by the too steep low-energy rise of the soft 
power-law. This soft power-law component, in addition, causes an infrared divergence after 
eliminating the photo-absorption. Further considering the unacceptable reduced chi-square 
value, this Model B is unlikely to represent the soft component. Although the two-blackbody 
model (Model C) can gives a better fit, the chi-squared value is still somewhat large. In 
addition, the associated hard-tail shape, Fh ~ 1.5, is considerably larger than that obtained 
with the HXD alone. 

In contrast to Model B and Model C which employ an ad-hoc combination of simple 
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model components, the two more physical models, Model D and Model E, were found to be more 
successful with reduced chi-square values of 1.08 and 0.88, respectively (table 1). Both models 
assume a scenario that some seed photons are repeatedly up-scattered by hot electrons in the 
stellar magnetosphere; the former is empirically defined (Enoto et al. 2010a), while the latter 
is based on a simplified ID semi-analytical modeling (Rea et al. 2008). Thus, these successful 
wide-band fits suggest that the known deviations of the soft component from a simple blackbody 
can be ascribed, at least partially, to up-scattering of soft seed photons off hot electrons located 
somewhere. Below, let us discuss a possible scenario, employing the latter modeling (Model 
E), which is more specific. 

As shown in fiture 6, Model E has given an acceptable fit to the 0.8-200 keV spectrum, 
with a seed-photon temperature of kT = 0.30, an optical depth of r = 2.0, and an electron 
thermal velocity of /3 = 0.3. These parameters are consistent with kT = 0.30, r = 1.9, and 
P = 0.33, derived by a combined fitting of the 2004 XMM-Newton observation and 2003-2006 
INTEGRAL observations using the same RCS plus a hard tail power-law modeling (Rea et 
al. 2008). Thus the soft X-ray spectral shape measured in 2007 has little difference from that 
in 2004. Therefore, the persistent soft X-ray emission of 4U 0142-1-61 is relatively stable at 
least on a time scale of a few years. In addition, the absorbing column density measured with 
Suzaku, A^H = 6.72l|5:o4 x 10^^ cm'^, is in a good agreement with A^h = (6.4 ±0.7) x 10^^ cm^ 
derived by absorption edges of elements such as O, Fe, Ne, Mg, and Si (Durant & van Kerkwijk 
2006b). 

When using Model E, the overall soft- component shape mainly determines the blackbody 
temperature kT, while the data excess above the blackbody toward higher energies, to be called 
"soft tail", specifies r and /3. The derived blackbody temperature, kT ~ 0.3 keV, is similar to 
those derived by Model C and Model D. (If using a more familiar Comptonized blackbody 
model, "compbb", we obtain r ~ 2.9, kT ~ 0.32 keV, and an electron temperature kT^ ~ 2.3 
keV.) Then, let us discuss whether these conditions of Model E are physically realistic in view 
of general pictures of magnetized neutron stars. These parameters imply that seed photons 
with kT ~ 0.3 keV, presumably from the stellar surface, pass through a scattering slab with an 
optical depth of ~ 2, filled with electrons which have typical velocities of 0.3c, where c is the 
speed of light. Since the magnetic field of 4U 0142+61 is estimated as -B ~ 1.3 x 10^^ G from 
the period P ~ 8.68 s and its derivative P ~ 2 x 10~^^ s s~^, the corresponding Goldreich- Julian 
density (Goldreich & Julian 1969) becomes ncj = 7 x 10i3(P/l sy\B/10^^ G) cm'^ = 1 x 10^^ 
cm~^. If assuming this density, the Thomson cross-section ctt, and a slab thickness of / ~ 10'' 
cm, the corresponding optical depth becomes r ~ ?igjC"t^ ~ 7 x 10~^, which is much smaller 
than the value, r ~ 2, required by the soft-tail feature in the data. However, consideration 
of more detailed physics of magnetars, as employed in the RCS model, allows us to solve 
this discrepancy by invoking two effects; a larger electron density and an enhanced electron- 
photon cross-section. The former is expected under a twisted magnetic configuration, where 
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electric currents are induced to give ~3-4 orders of magnitude higher electron density than 
ncj (Thompson et al. 2002). The latter, an enhanced cross-section, can occur at a cyclotron 
resonance energy in the scattering slab (Lyutikov & Gavriil 2006). These larger density and 
cross-section are consistent with a basic scheme of the magnetar hypothesis. 

4-3. Hard X-ray emission of 4U 0142+61 

The fit with Model E has given the photon index of the hard component as Fh = 
0.89tQ}Q(stat.)to lo(sys-)- This is consistent with, and of similar quality to, the value of 
0.93 ±0.06 measured by INTEGRAL (den Hartog et al. 2008), or that based on the XMM- 
Newton and INTEGRAL data, 1.1 ±0.1 (Rea et al. 2008), even though the Suzaku exposure, 
~100 ks, is much shorter than that with INTEGRAL. Among the magnetars shown in figure 4 
in Enoto et al. (2010c), 4U 0142±61 thus exhibits one of the hardest power-law components. 

The emission mechanism of the hard X-ray component is a topic of extensive discussion. 
Although several scenarios have been proposed, none of them succeeds to explain all the obser- 
vational properties. So far, one of the biggest difficulties therein has been how to explain the 
extremely hard photon indices (Fh ~ 1 or less; figure 7). Now, the difficulty is enhanced by the 
systematic change in Fh (from ~ 1.7 to ~0.4) in correlation with the characteristic age (Kaspi & 
Boydstun 2010; Enoto et al. 2010c). Indeed, this effect cannot be readily explained by the non- 
thermal radiation model from fast-mode break down (Heyl & Hernquist 2005), or the thermal 
bremsstrahlung model from a transition layer (Thompson & Beloborodov 2005; Beloborodov 
& Thompson 2007), or the resonant inverse-Compton up-scattering model (Baring & Harding 
2007). 

A possible observational hint to the emission mechanism of the hard component may be 
provided by the highest energies to which the hard tail extends. There is a relatively strong 
COMPTEL upper limit above ^750 keV (Kuiper et al. 2006). In order to make this point 
clear, we show, in figure 8, a multi-band spectral energy distribution of 4U 0142±61, incorpo- 
rating previously reported multi- wavelength observations. Thus, if time variations are negligible 
among these observations, the hard X-ray spectrum must turn over at energies between ~200 
and ~750 keV. The Suzaku lower limit on the cutoff, Ecnt > 180 keV, is consistent with the 
above estimate. If this hard power-law component originates from non-thermal emission from 
accelerated particles, the cutoff would reflect the maximum particle energy. However, consid- 
ering the difficulty with the ordinary non-thermal emission processes by energetic particles to 
explain the values of Fh and its systematic source-to-source differences, the cutoff may alter- 
natively be interpreted as an effect of the extremely high magnetic field, typically exceeding 
the critical field B^r = 4.41 x 10^'^ G, where gamma-rays cannot directly escape due to photon 
splitting (Harding & Lai 2006). 

This photon-splitting view implies a possibility of the hard X-ray emission arising via 
down-cascade of energetic gamma-rays (Harding et al. 1997; Enoto et al. 2010d). Near the 
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stellar surface, gamma-rays above a few hundred keV can be produced through, e.g., resonant 
cyclotron scattering, annihilation of electron-positron pairs at the dense stelar atmosphere, and 
sporadic reconnection. In ultra-strong magnetic fields exceeding Bcr, electron-positron pair 
cascades are suppressed (Baring & Harding 2001), while the photon sphtting may be dominant 
(Baring & Harding 2001; Harding & Lai 2006). As a result, these gamma-rays from the surface 
may repeatedly split into lower energy photons, since the photon splitting process has no low- 
energy threshold unlike the one-photon pair production process. In fact, Harding et al. (1997) 
reported that the photon splitting becomes important at 5 ^ O.SBcr, since the attenuation 
length due to splitting become comparable to or less than those for pair production . This 
process can also explain the differences in Fh among magnetars, in such a way that higher fields 
of younger objects will allow the photon-splitting cascade to proceed down to lower energies, 
and hence to make the continuum softer. 

4.4- Energetics 

Assuming an isotropic emission at an distance of 3.6 kpc (Durant & van Kerkwijk 2006a), 
the soft (RCS) and hard components are implied to have unabsorbed luminosities of 2.8 x 10^^ 
erg s~^ and 6.8 x 10'^^ erg s"^, in the 1-10 keV and 10-70 keV band, respectively. If calculated in 
the 2-10 keV band, our soft-component luminosity becomes 1.1 x 10'^^ erg s~^. Since this agrees 
with the previous measurements (Durant & van Kerkwijk 2006a), the soft component is inferred 
to be stable, making some theoretical prediction (Thompson & Duncan 1996) consistent with 
the observation. Adding up the soft and hard components, the absorption-corrected 1-200 keV 
luminosity becomes 5.2 x 10"^^ erg s~^, which can be decomposed into ~54% and ~46% carried 
by the soft and hard components, respectively. Since the spin-down luminosity of the source is 
only 1.2 X 10^^ erg s~^, neither components can be powered by the rotational energy. 

4-5. Pulse profiles 

The soft and hard X-ray pulse profiles, obtained with Suzaku (figure 4), are generally 
consistent with those measured previously with INTEGRAL and RXTE (Kuiper et al. 2006; den 
Hartog et al. 2008). In energies below a few keV, the profile has two peaks, and one of them 
(phase at ~0.5 in our figure 4 and 0.1 in figure 5 of Kuiper et al. 2006) becomes weaker towards 
10 keV, but it partially recovers towards a few tens keV at a somewhat smaller pulse phase. 

As seen above, the pulse profiles of 4U 0142-1-61 depend in a complex way on the energy. 
However, the main peaks of the two spectral components are at the same phase ~1 in figure 4, 
and their emission regions seem to be located at the same rotation phase. Regarding the soft 
component as an optically-thick thermal emission from the polar-cap regions of the neutron 
star, we further speculate that the spectral hard-tail component, which are generally pulsed 
strongly (e.g., Kuiper et al. 2006), is also emitted from the polar-cap regions. This inference 
gives an additional support to the photon-splitting mechanism, because the input gamma-rays 
will be produced, e.g., via electron-positron annihilation, mainly at the polar-cap regions. 
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In subsection 3.4 and figure 5, we presented evidence for sliort-term (in 30 ks) variations 
in tlie pulse profile (or overall signal intensity). During the last 30 ks of our exposure, when 
the hard X-ray intensity slightly decreased (figure 1), the HXD-PIN pulse profile possibly 
became more strongly pulsed and more double peaked. These effects, if real, may reflect 
some sporadic processes involved in the hard-tail production mechanism. For example, if the 
persistent emission of a magnetar is formed by an assembly of numerous small short bursts 
(Nakagawa et al. 2007), the persistent intensity should fluctuate due to statistical fluctuations 
of the number of small bursts. We expect the soft spectral component to be more stable, 
because of the heat capacity of the neutron star. More quantitative evaluation of these issues 
will be presented elsewhere. 

5. Conclusion 

We observed the prototypical anomalous X-ray pulsar 4U 0142+61 with Suzaku, and 
quantified its wide-band spectra, spanning from 0.4 to ~70 keV or higher (subsection 3.2). 
This is a first detailed broad-band study of this source performed using a single satellite within 
a short exposure (~100 ks), without cross-calibration uncertainties among multiple satellites. 
We reconfirmed the hard-tail component above 10 keV, with high significance of ^ 3cr. Detailed 
analyses of the soft and hard X-ray spectra, simultaneously fitted by several spectral models 
(subsection 3.5.1), shows that the soft component is represented by Comptonized blackbody 
model or Resonant Cyclotron Scattering model with a seed-photon temperature of kT ~ 0.3 
keV. The hard-tail component can be fitted by a power-law with a photon index of Fh ~ 0.9, 
together with a constraint on its cutoff as E^ut > 180 keV. The folded pulse profiles in the 
12-40 keV range exhibited evidence for shape or intensity variations on a time scale of 30 ks 
(subsection 3.4), suggesting some sporadic process in the hard X-ray production process. We 
propose the photon splitting process as a promising production mechanism for the hard X-ray 
component. 

We are grateful to all the Suzaku team members, including in particular M. Morii for 
his elucidating discussion. 
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Table 1. Spectral parameters of 4U 0142+61 observed in 2007 August*. 
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* BB, PL, CBB and RCS represent blackbody, powe-law, (magnetic) comptonized blackbody, and resonant 
cyclotron scattering, respectively. All the quoted errors arc only statistical at the 90% confidence level. 
Lower-case characters "h" and "s" in model names represent "hard" and "soft", respectively, 
t The absorbed X-ray fluxes (IQ— ergs s~^ cm"^ ) in the 2-10 keV. 
t The X-ray flux (lO'^^ ergs s^^ cm"^ ) tj^e 15-60 keV. 
§ Distance to 4U 0142+61 is assumed at 3.6 kpc. 

I The parameter a is related with the soft tail power-law as a = F; 
process. Details arc given in an appendix of Enoto et al. (2010a). 
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Fig. 1. (a) Background-subtracted 0.4-10 keV light curve of XIS3, witli a binning of 160 s. (b) The 
NXB-subtracted HXD-PIN light curve in the 10-70 keV range, with 550 s binning. It still includes the 
CXB contribution by ^0.03 c s~^, which is shown as a dotted line. 
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Fig. 2. (a) Background-subtracted Suzaku spectra of 4U 0142+61 obtained with the XIS, HXD-PIN, 
and HXD-GSO. Not only the NXB but also the CXB contribution was subtracted from the HXD-PIN 
data. The error bars are statistical only for the XIS and HXD-PIN, while those of HXD-GSO include 
systematic errors as well, (b) The same spectra as presented in panel (a), but divided by those of the 
Crab Nebula. 
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Fig. 3. Pcriodograms from XIS0+XIS1+XIS3 (panel a: 0.8-10 keV) and HXD-PIN (panel b: 12-40 
keV), calculated using 0.54 s binned light curves and 16 bins for one pulse period. The backgrounds are 
inclusive. 



19 



O 0.58 



. 1 1 III 1 ■ 
■ 1 . L/t^ . . 1 


■ rtj 1 ' ' ' ' 1 ' ' ' rtj ■ 
. . 1 . L/t^ . . 1 . . . 


: \P, \ 


' 







-i . . . . 1 . 


. . 1 .... 1 





(a) 



(b) 



(c) 



(d) 



1 

Phase 



Fig. 4. Background-inclusive XIS (XISO, XISl, and XIS3 summed), HXD-PIN, and HXD-GSO pulse 
profiles, obtained by folding the total data at the pulse period of equation (1), in (a) 0.8-4.0, (b) 4.0-10, 
(c) 10-70, and (d) 80-150 keV ranges. The HXD data is corrected for the dead time. Corresponding 
background rates are estimated to be (1.34 ±0.04) x 10~^ counts scc~^, (1.94 ±0.14) x 10~^ counts sec""'^, 
(5.49±0.11) X 10^^ counts sec^^, and 5.00±0.21 counts sec~^, for panels (a), (b), (c), and (d), respectively. 
These background rates include only statistical errors for the XIS and both statistical and systematic 
errors for the HXD. We assumed 2% and 0.6% systematic errors for the PIN-NXB and the GSO-NXB, 
respectively. 
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Phase 



Fig. 5. NXB- and CXB- subtracted and dcadtimc-corrccted 12-40 kcV HXD-PIN pulse profiles, during 
(a) 1st, (b) 2nd, and (c) 3rd 30 ks intervals, respectively. Error bars are statistical (la) only, whereas the 
sytcmatic errors of the background subtraction is ^0.004 c s~^. 
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Fig. 6. (a) The same Suzaku spectra of 4U 0142+61 as shown in figure 2, fitted by RCS (green)+PL 
(magenta) model, namely Model E. (b) Residuals from the fit in panel (a). (c)(d)(e)(f) The same as panel 
(b), but the soft component is represented by BB (Model A), BB+PL (Model B), 2BB (Model C), and 
CBB (Model D) models, respectively. 
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Fig. 8. A multi-band spectral energy distribution of 4U 0142+61. The Suzaku data are shown in the 
0.8-10 keV and 12-200 keV energy ranges from the XIS and the HXD, respectively. The soft and hard 
components of Model E are represented by dashed lines after corrections for the interstellar absorption, 
and are extrapolated to lower energies. The previous X-ray results are also shown; the 0.1-200 keV band 
from the XMM-Newton/PN (Rea et al. 2007; Abdo et al. 2010) and the INTEGRAL/ISGRI data sets 
(den Hartog et al. 2008). In the gamma-ray energy band, the 2a CGRO/COMPTEL upper limits (den 
Hartog et al. 2006; Kuiper et al. 2006) and 95% Fermi/LAT upper limits (Abdo et al. 2010) are ploted. 
In the infrared and optical ranges, Spitzer observations at 4.5 /im and 8.0 fim (Wang et al. 2006) and 
some data sets from Gemini and Keck II (Durant & van Kerkwijk 2006c; HuUeman et al. 2004) are shown, 
where the optical data are de-extincted assuming a reddening value of Ay ~ 3.5. In the radio band, a 2cr 
upper limit from 1.38 GHz WSRT continuum (den Hartog et al. 2007), and the 1.4 GHz VLA upper limit 
(Gaensler et al. 2001) arc shown. 
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